Docket No. 210731US2 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

IN RE APPLICATION OF: Jean- Yves DELABBA YE, et al. GAU: Q ^= 

SERIAL NO: New Application EXAMINER: £ , W= 

FILED: Herewith c 

FOR: DEVICE FOR THE ANALYSIS OF ELECTROMAGNETIC SIGNALS =>|p? 

REQUEST FOR PRIORITY gSTS^ 

ASSISTANT COMMISSIONER FOR PATENTS ^ == 

WASHINGTON, D.C. 20231 

SIR: 

□ Full benefit of the filing date of U.S. Application Serial Number, filed, is claimed pursuant to the provisions of 
35 U.S.C. §120. 

□ Full benefit of the filing date of U.S. Provisional Application Serial Number , filed , is claimed pursuant to the provisions of 
35 U.S.C. § 119(e). 

■ Applicants claim any right to priority from any earlier filed applications to which they may be entitled pursuant to the 
provisions of 35 U.S.C. §119, as noted below. 

In the matter of the above-identified application for patent, notice is hereby given that the applicants claim as priority: 

COUNTRY APPLICATION NUMBER MONTH/DAY/YEAR 

FRANCE 00 09 1 3 1 July 1 2, 2000 



Certified copies of the corresponding Convention Application(s) 
■ are submitted herewith 

□ will be submitted prior to payment of the Final Fee 

□ were filed in prior application Serial No. filed 

□ were submitted to the International Bureau in PCT Application Number . 

Receipt of the certified copies by the International Bureau in a timely manner under PCT Rule 17.1(a) has been 
acknowledged as evidenced by the attached PCT/IB/304. 

□ (A) Application Serial No.(s) were filed in prior application Serial No. filed ; and 
(B) Application Serial No.(s) 

□ are submitted herewith 

□ will be submitted prior to payment of the Final Fee 




22850 

Tel. (703)413-3000 
Fax. (703)413-2220 
(OSMMN 10/98) 



Respectfully Submitted, 



OBLON, SPIVAK, McCLELLAND, 
MAIER & NEUSTADT, P.C. 

Marvin J. Spivak 
Registration No. 24,9 1 3 

C. Irvin McClelland 
Registration Number 21,124 



This Page Blank (uspto) 



IN 5 ! 



I INSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



Of 



o 



BREVET D'lN VE N T I O N 



ti 



CERTIFlS^f D'UTILITE - CERTIFICAT D f Kl%lflON 



CQPlfe OFPj(D§LLE 




, certifite4#fo1r» prop 



Fait a Paris, le 



\ 5 JU1« 2001^ 



i \ 

: \ 

u 
j! 

s 

// 



Pour le Directeur general de Hnstitut' 
national de la propriete industrielle 
Le Chef du Departement des brevets 




Marline PLANCHE 



INSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



SIEGE 

26 bis, rue de Saint Petersbourg 
75800 PARIS cedex 08 
Telephone : 01 53 04 53 04 
Telecopie : 01 42 93 59 30 
http://www.tnpi.fr 



THIS PAGE BLANK cuspio) 




NATIONAL Bt 

Uf MMill l 

l«»0tTKWU.B 

26 bis. rue de Saint Petersbourg 
75800 Paris Cedex 08 

Tel6phone : 01 53 04 53 04 Telecopie : 01 42 94 86 54 



BREVET D INVENTION 
CERTIFICAT DUTILITE 

Code de la propriete intellectuelle - Uvre VI 

REQUltTE EN DtUVRMlCE 1/2 



N° 11354*01 



^T EDES ^JUIL 2000 
75 IN PI PARIS 



DATE 
UEU 



N° DENREGJSTREMENT 
NATIONAL ATTWBUE PAR UINPI 
DATE DE DEPOT ATTRIBUEE 
PAR L'INPI 



Confirmation orun aepoi pa^eiwwpw 
Q NATURE DE LA DEMANDE 


Cochez Tune des 4 cases suiwantes | 


Demande de brevet 


□ 1 


Demande de certificat d'utilite 


□ 1 


Demande divisionnaire 


□ 








Demande de btvtvt initiale 




Date 


/ 


/ 1 


ou demande de certifwat d'utilite initiale 


N° 


Date 


/ 


' ' 1 


Transformation d'une demande de 
brevet europeen Demande de brevet initiale 


□ 

N° 


Date ! 


/ 


/ 1 



- 4 R6serv6 k HNPI \ 



0009131 



Vos references pour ce dossier 

(facultatif) THOMSON Aff. 594 (16459) 



Cet impri me est a remplir lisiblement a I'encre noire 

H 



08 540 W /26G899 



NOM ET ADRESSE DU DEMANDEUR OU DU MANDATAIRE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 



CABINET NBTTER 
40 rue Vignon 
75009 PARIS 



XITRE DE LMNVENTION (200 caracteres ou ©spaces maximum) 

Dispositif pour Tanalyse de signaux electromagnetiques. 



Q DECLARATION DE PRIORITE 
OU REQUtTE DU BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANT^RIEURE FRANQAISE 


Pays ou organisation 1 
Date / / N° I 

Pays ou organisation 1 
Date / / N° 1 

Pays ou organisation I 
Date / / N° 

□ S'il y a d'autres prlorites, cochez la case et utilisez rimprime aStritea | 


Q DEMANDEUR 


□ S'il y a d'autres demandeurs, cochez la case et utilisez rimprime aSuiten 1 


Nom ou denomination sociale 


THOMSON-CSF 1 


Prenoms 




Forme juridique 


Societe Anonyme 1 


N° SIREN 




Code APE-NAF 


i '. '■ '■ ■ ■ ■ 1 


Adresse 


Rue 


173 boulevard Haussmann I 


Code postal et ville 


75008 1 PARIS 1 


Pays 




France 1 


Nationality 


francaise 1 


N° de telephone (facultatif) 




N° de telecopie (facultatif) 




Adresse electronique (facultatif} 






mOTITUT 
LO POOPWOTtt 

ineuoiaiati-a 



o 

CERfFBIFDCOT DWOOTlfe 
REQUITE EN D&JVRANCE 2/2 



H Reserve ai'lNPlf 



REMISE DES PIECES ____ 

date U2 JUDL 2000 
75 INPI PARIS 



LIEU 



N° D'ENREGISTREMENT 
NATIONAL ATTWBUE PAR LINPI 



©0©S>131 



Dfl 540 W /260899 



^©s tpAffttrencos pour ce dossier : 
(faatttatif I 



Norn 



Prenom 



Cabinet ou Societe 



N °de pouvoir permanent et/ou 
de lien contractuel 



Adresse 



Rue 



Code postal et ville 



N° de telephone tftiatltatif) 



N° de telecopie ifaadtatij) 



Adresse electronique tjaadtafij) 



THOMSON Aff. 594 (16459) 



PLA^AIS 



Jean-Yves 



Cabinet NETTER 



40 rue Vignon 



75009 | PARIS 



01 47 42 02 23 



01 47 42 60 02 



It] bgweotemr (sj 



Les inventeurs sont les demandeurs 



□ Oui 



E Non Pans ce cas ffoumroir mine desilgiraattioira dNmrontemts) separee 



RAPPQKT DE KECHECSCHE 



MmiSqiueinnteini* pour uoie demand© de forevett (y compirto dfoosfion et torcsfqirmatkw) 



Etablissement immediat 
ou etablissement differe 



Paiement echelonne de la redevance 



PaSememit an trofis veirsements, unlquMieiit pour les personnes physiques 

□ Oui 

□ Non „ — 



r£dujothoin dqj taux 
des kepevawces 



UJiiiqueinnent pour les personnes physiques 

□ Requise pour ta premiere fois pour cette invention ( jainchv tm arLuUnmiHrnfmitrntt ) 

□ Requise anterieurement a ce depot (jattuhv urn* cofiie de la decismu d'cutmisswit 
pom- cette imvntiait ott btdiquersa lyfeivttce) : 



Si vous avez utilise IMmipriime aSuifeen, 
indiquez ie noimbre de pages jointes 



tg SQGEWORE OU DESWAWDEUR 
00 OU RflAWDATADRE 
(Nom et qualfte du signataire) 



1197 (B) (M) 
QA1S 




VQSA DE LA PREFECTURE 
€>W OE UOWPO 




La loi n°78-17 du 6 janvier 1978 relative a I'informatique. aux fichiers et aux liberies s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de I'lNPI. 



IN 5 ! 



■ IKMITUT 
NATIONAL D( 
LA MOMim 
INOUtTOIILLI 



DtPARTEMENT DES BREVETS 

26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 01 53 04 53 04 Telecopie : 01 42 94 86 54 



BREVET D'INVENTION 
CERTIFICAT D'UTILITE 

Code de la propriete intellectuelle - Livre VI 

DESIGNATION D'INVENTEUR(S) Page N° . } . / . - 1 . 
(Si le demandeur n'est pas IMnventeur ou I'unique inventeur) 



N° 11235*02 



OB 113 W/260899 



Vos references pour ce dossier 

(facuitatif) 


THOMSON-CSF Aff . 594 


N° D'ENREGISTREMENT NATIONAL 


N° 00 09131 du 12 juUlet 2000 



TlTRE DE L'INVENTION (200 caracteres ou espaces maximum) 

Dispositif pour l'analyse de signaux electromagnetiques. 



LE(S) DEMANDEUR(S) : 

THOMSON-CSF 



DESIGNE(NT) EN TANT QU'INVENTEUR(S) : (Indiquez en haut & drolte «Page N° 1/1» S'il y a plus de trols inventeurs, 
utilisez un formulaire identique et numerotez chaque page en indiquant le nombre total de pages). 



Nom 


DBLABBAYB 


Prenoms 


Jean-Yves 


Adresse 


Rue 


13 avenue du President Salvador Allende 


Code postal et ville 


94117 ARCUBIL 


Societe d'appartei 


lance (facuitatif) 




Nom 


COUDERETTB 


Prenoms 


Andre 


Adresse 


Rue 


13, avenue du President Salvador Allende 


Code postal et vilte 


94117 ARCUBIL 


Societe d' a p parte 


nance (facuitatif) 




Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 




Societe d'apparte 


nance (facuitatif) 




DATE ET SIGNATURE(S) 
DU (DES) DEMANDEUR(S) 
OU DU MANDATAIRE 
(Nom et qualrte du signataire) 


Paris, le 18 mai 2001 

N° Conseil 92-1197 (B) (M) 



La loi n°78-17 du 6 janvier 1978 relative a I'informatique, aux fichiers et aux liberies s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de I'INPI. 



DE594 „ FRD 



Dispositif pour 1 1 analyse de sicrnaux 61ectromagn 6tiques<, 

L 1 invention concerne 1 1 analyse de signaux 61ectromagn6tiques 
a priori inconnus • 

De fagon g6n6rale, on sait convertir un signal 61ectromagn6- 
tique regu ("signal r6el" ) en un signal complexe, repr^senta- 
tif en amplitude et phase du signal 61ectromagn6tique regu, 
Ceci se fait en pratique apr&s avoir r6duit la frequence, pour 
descendre g6n§ralement jusqu'en bande de base* Le signal 
complexe poss&de deux composantes, l'une en phase et 1' autre 
en quadrature, h la maniere des nombres complexes, d'ou son 
nom. 

Les techniques actuellement disponibles pour 1' analyse de 
signaux inconnus consistent Si rechercher dans le signal regu 
des attributs invariants, par translation de frequence ou de 
niveau, notammento Elles peuvent faire intervenir une 
transformation de Fourier, ou bien une transformation de type 
"ondelettes" , ou encore une transformation de WIGNER-VILLE , 
par exemple* 

Bien que l'^tat de la technique soit difficile & percevoir 
exactement en la matiSre, la Demanderesse considSre actuelle- 
ment comme connu de rechercher une estim6e de la moyenne de 
la "frequence instantan^e" , et d'effectuer une demodulation 
du signal par cette frequence o Cela pourra etre utilise dans 
la mise en oeuvre de 1' invent ion „ 

La Demanderesse s'est pos# le probl£me d'obtenir des moyens 
qui aident k trouver le message contenu dans un signal, ou 
meme permettent de trouver compl&tement ce message, lorsqu'on 
ne sait rien sur la modulation exacte et 1 ■ alphabet de 
symboles utilises h 1' Emission, quand, par exemple, la 
modulation est num6rique« 

Elle s'est int6ress6e & certains types g€n6raux de modulation, 
en particulier & la modulation num^rique linSaire, pour 



lesquels elle a trouv<5 des solutions, qui sont les elements 
constitutifs de la pr6sente invention. 



L ' invention porte done, notamment, sur un dispositif d'aide 
a 1' analyse de signaux contenant une modulation num6rique par 
des symboles. Ce dispositif comprend: 

- une memoire de signal, pour stocker un signal numerique 
complexe <z(t)), repr^sentatif en amplitude et phase d'un 
signal capt£, sur une dur6e choisie, et 

- des moyens de traitement, agenc^s pour rechercher dans le 
signal complexe des propriStes relatives a sa frequence 
porteuse et a la modulation de celle-ci, en fonction d'un 
modele de modulation, en particulier linSaire. 

Selon un premier aspect de 1' invention, les moyens de 
traitement comprennent des moyens pour determiner une estimee 
du rythme (1/T) de la modulation, et des moyens de projec- 
tion, agenc<§s pour calculer les composantes (z p (k)) du signal 
complexe dans une base de fonctions (^(t)), laquelle est 
param6tree selon ledit rythme (1/T) de la modulation, et des 
moyens de calcul, operant sur ces composantes, afin de 
determiner au moins une estimee relative a au moins une 
propriety du signal complexe, dans le groupe de propri<§tes 
comprenant la forme d' impulsion eiementaire (g(t)) du signal 
complexe, la suite des symboles (a(k)) du signal complexe et 
la porteuse f 0 utilised. 

Avantageusement : 

- la base de fonctions est paramStree pour presenter au moins 
deux echantillons par periode (T) de la modulation, et plus 
generalement q Echantillons par periode (q e Net q * 2); 

- elle comprend au moins deux fonctions (^(t) , /*> b ( t >) se 

d^duisant l'une del' autre par une translation temporelle de 
periode choisie (T/2 ou plus generalement T/q, q e N et q*2); 

- elle peut §tre fondee notamment sur des fonctions rectangu- 
laires tempore llement adjacentes les unes aux autres, ou 
encore sur des fonctions du type "cosinus sur61ev6" . 
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Dans un mode de realisation actuellement prefere. les moyens 
de projection comprennent: 

_ des moyens definissant un filtre numerique, ayant une 
reponse impulsionnelle sensiblement egale a 1'une desdites 
fonctions, ce filtre numerique recevant le signal complexe, 

f des moyens d'echantillonnage numerique repete de la sortie 
de ce filtre a une cadence choisie (2/T ? q/T, q e N et q £ 2 ) , 
suivant le type de base choisie. 

Selon un autre aspect de 1' invention, les moyens de traite " 
ment sont agences pour determiner une estimee approchee f a de 
la frequence porteuse f 0 du signal complexe, ainsi que pour 
demoduler ce signal complexe par cette estimee f a . Les moyens 
de projection sont alors agents pour operer sur le signal 
complexe (z(t)) apres sa demodulation par cette 
approchee, tandis que ladite base de fonctions est de basse 
frequence, sensiblement comme le spectre du signal complexe 
qui, apres demodulation, est lui aussi basse frequence. 

Selon encore un autre aspect de 1'invention, les moyens de 
calcul comprennent des moyens de calcul matrxcxel sur 
lesdites composanteso 

De preference, on calcule d'abord une solution Particuli^re 
du couple d'inconnues (g p (t), a(k) ) du modele du sxgnal apres 
projection, sous la forme d'une fonction de support mxnxmal 
(h (t)), assortie d'un train de symboles (c m (k) ) . De ce 
r^idu, on pourra tirer un residu de filtrage (a,) *ux, 
associe a la forme d' impulsion de support mxnxmal, donnera 
une estimee de g(t). On verra qu'on peut en tirer ensuxte 
symboles (a(k)) du signal initial. 

Exprimee ci-dessus en termes de dispositif , 1'invention peut 
aussi §tre definie sous la forme de precedes. 

D'autres caracteristiques et avantages de 1'invention 
apparaitront a 1'examen de la description detaillee qux va 
suivre, ainsi que des dessins annexes, sur lesquels : 



- la figure 1 illustre schematiquement , sous forme de 
diagramme-blocs , le principe du passage en signal complexe 
(100); 

- la figure 1A illustre schematiquement un premier mode de 
realisation permettant de mettre en oeuvre le principe du 
passage en signal complexe de la figure 1; 

- la figure IB illustre schematiquement un second mode de 
realisation permettant de mettre en oeuvre le principe du 
passage en signal complexe de la figure 1; 

- la figure 2A illustre un exemple de fonction de modulation 
elementaire sous la forme d'une impulsion en "racine de 
cosinus surelev6'\ en fonction du temps expriiruS en periodes 
T; 

- la figure 2B illustre une superposition de signaux a(k)g(t- 
kT) (ou symboles) emis success ivement , en fonction du temps 
exprim6 en periodes T; 

- la figure 2C illustre la resultante des signaux de la 
figure 2B; 

- la figure 2D illustre une repartition de symboles (ou 
constellation) dans le plan complexe, dite "alphabet"; 

- la figure 3 illustre les principales etapes (ou stades) de 
traitement du procede selon 1' invention; 

- les figures 4A a 4D illustrent respectivement une reponse 
impulsionnelle d'un filtre passe-bas reel classique, une 
fonction de transfert du filtre de la figure 4A, un spectre 
de signal d' entree, et le spectre du signal complexe du 
filtre ayant la fonction de transfert de la figure 4B; 

- les figures 5A a 5D illustrent respectivement le spectre 
d'un signal obtenu a partir d'une FFT, un filtrage non 
lineaire de la sortie de la FFT de la figure 5A, un exemple 
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de traitement des raies de la figure 5B, et un exemple 
d' estimation de la raie en + 1/T de la figure 5C; 

- la figure 6 illustre sous forme d'etape un exemple de 
traitement permettant d'obtenir la valeur 1/T de la raxe 
isolee de la figure 5D; 

- les figures 7A a 7D illustrent respectivement un exemple de 
signal defini par un nuage de points dans la plan complexe, 
le spectre du signal complexe de la figure 7A, le spectre du 
signal complexe de la figure 7B apres demodulation, et la 
densite spectrale du signal de la figure 7A avant et apres 
demodulation; 

- les figures 8A et 8B illustrent respectivement les projec- 
tions de h. et du signal dans une base de fonctions choxsxe; 

- la figure 8C illustre une approximation des fonctions g B (t) 
dans la base de fonctions choisie aux figures 8A et 8B; 
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figure 9 illustre sous forme d'etapes un exemple du 
deroulement des traitements visant a estimer h m (t); 

- les figures 10A a 10C illustrent respectivement le spectre 
de la fonction de transfert du filtra de reponse impulsxon- 
nalla pour un signal centre du type de celui de la figure 7C 
et pour un premier type de base, le spectre de la fonction * 
transfert de fonctions <p en cosinus sureleves, et le spectre 
de la fonction de transfert du filtra de reponse impulsxon- 
nelle correspondant a la figure 10B; 

la figure 11 illustre 1 ' approximation de h^t) par sa 
projection h pm (t) dans la base de fonctions de la figure 10C; 

la figure 12 illustre la repartition des valeurs propres de 
la matrice M L1 , rangees par ordre decroissant de grandeur en 
abscisse, lorsque la projection correspond a un fxltrage 
suivi d'un echantillonnage a q = 2 points par symboles; 



_ la figure 13 illustre la repartition des valeurs propres de 
la matrice M L1 , rangees par ordre decroissant de grandeur en 
abscisse, lorsque la projection correspond a un filtrage 
suivi d'un echantillonnage a q > 2 points par symboles; 

- les figures 14A a 14C illustrent respectivement , en 
fonction du temps exprime en periodes T, la fonction h pm (t), 
le residu de filtrage, et la fonction elementaire de modula- 
tion modul^e g pm (t); 

- la figure 15 illustre sous forme d'etapes un exemple du 
deroulement des traitements visant a estimer les symboles 
a^k) a partir des trains de symboles possibles; 

- la figure 16 illustre un exemple de dispersion du fait du 
bruit des symboles a(k), dans le plan complexe et pour une 
plage de frequence (residu de frequence) bien centree; 

- la figure 17 illustre un second exemple de dispersion du 
fait du bruit des symboles a(k), dans le plan complexe et 
pour une plage de frequence mal centree; 

- la figure 18 est une representation tri-dimensionnelle de 
la densite de probability des symboles a(k) de la figure 16; 

- la figure 19 est une representation tri-dimensionnelle de 
la densite de probability des symboles a(k) de la figure 17; 

- la figure 20 illustre sous forme d'etapes un exemple du 
deroulement des traitements visant a estimer la position des 
ytats possibles; et 

- la figure 21 illustre sous forme d'etapes un exemple du 
deroulement des traitements visant a estimer la suite des 
symboles ef f ectivement emise, a partir du train de symboles 
bruites a(k) et des etats possibles des symboles e ± calcules 
en reference a la figure 20. 
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L es dessins annexe so„ t . pour dessentiel - 
certain. En conseguence, ils pourront non seulement server > 
mieux faire comprendre la description, maxs aussx contrxbuer 
a. la definition de 1' invention, le oaa echeant. 

us expressions Mathematics elaborees, ou gui 
repetitxvement, sont regroupees dans une annexe 
partie de la description. Elles sont designees par la lettre 
E, suivie d'un nombre de reference. 

Oans ces expressions, le vecteur unitaire des imaginaires 
complexes est classxquement note j, avec 3 - - • 
ailleurs, le symbole -i», mis en indice ou comme variable de 
section, est utilise * plusieurs 

significations differentes; en effet, le nombre important de 
relation, invokes ne permet pas d^tiliser 

une lettre differente * chaque fois, et le contexte leve en 
princIpe toute ambiguxte. De meme, les ensembles des valeurs 
prises par les variables x, k. t.... ne sont pas generalement 
Lplicites, sauf lorsgue gue cela s-avere n.cessa.re Ces 
ensembles se deduisent du contexte en prenant toutes les 
vale-a possibles sur la totalite de la duree du sxgna 
analyse. L es expressions des valeurs minimales - —-ales 
en fonction des inconnues des modeles ne presente pas 
"interet pratigue dans l.algoritnme. En outre, ^ 
en exposant indigue le conjugue d-un nombre complexe s xl 
Tzll d-une matrice, il indigue sa transposee con 3 ugu*e. 
Enfin, deux caracteres adjacents peuvent indiguer un product 
si le contexte indigue gu-ils relevent de deux arables 
differentes = par exemple, le produit de 2 par 11 est note 
2L1 ou 2oLlo 

L es mots "estimer- et ..estimation" visent devaluation d- une 
grandeur inconnue, sans aucun a priori sur la precxsxon de 

35 cette Evaluation. 

Enfin, pour alleger 1'ecriture et la description, le bruit 
I est pas pris en compte dans les expressions mathematxgues 
donn^es ci-apres. 
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L' invention concerne 1' analyse de signaux de communication. 
Dans ce domaine, on part d'un signal de communication capte, 
dont on s'ef force de determiner le contenu informatif brut 
(c'est-a-dire independamment d'un eventuel chif f rement ) . On 
suppose que le signal capt6 a 6t<5 enregistrS, et qu'il est 
analyse" en temps diff6r6. 

Selon l'ouvrage "Digital Communications", John G. PROAKIS, 
(Chapitre "Representation of digitally modulated signals", p. 
163, seconde Edition, Mc Graw Hill Book Company), un signal 
correspondant a une modulation numerique lineaire se mod«§lise 
de la maniere indiquee dans la relation El (voir 1' annexe ou 
sont regroupees la plupart des formules), relation dans 
laquelle : 

- les a(k) (k entier positif ) sont les syroboles d' information 
transmis, qui ne peuvent prendre qu'un nombre fini de valeurs 
possibles, ou 6tats, 

- g(t) est la fonction de modulation 61ementaire (ou "forme 
d' impulsion dlementaire" ) , appliquee a chaque symbole, 

- 1/T est le rythme d' emission des symboles, et 

- f 0 est la frequence porteuse utilised . 

La fonction de modulation Slementaire (c'est-a-dire appliquee 
pour chaque symbole) est par exemple une forme d' impulsion en 
"racine de cosinus sur«§leve- (Fig. 2A) . Une telle forme 
d'onde est decrite par exemple dans l'ouvrage precitd de 
Proakis, page 536. Le plus souvent, les fonctions de modula- 
tion Slementaire rencontrees en pratique, ont, comme dans 
1' exemple ci-dessus, des longueurs plus grandes que T. 11 en 
decoule que les differents signaux a(k)g(t-kT) (ou symboles) 
emis successivement toutes les T secondes se superposent 
partiellement (Fig. 2B). D'apres la resultante de cette 
superposition (Fig. 2C), il est clair que ce phenomene 
contribue largement a rendre difficile 1' analyse du type de 
modulation. 

Par ailleurs, les symboles utilises sont choisis dans un 
alphabet ou constellation de symboles (exemple en figure 2D), 
eux aussi exprim^s en notation complexe. 



La Demanderesse s'est pose le probleme de determiner quel est 
le type de modulation employe dans un message capte, pour 
determiner ensuite le. parametres utilises pour cette modula- 
tion, jusqu'a 1'etat (ou valeur) de la suite de symboles 

5 ef f ectivement transmiSo 

La presente invention vient proposer une solution applicable 
notamment au modele de modulation numerique lineaire Svoque 
ci-dessus, lequel sera utilise dans les exemples donnas cx- 
LO apres. II. est rappele que font partie notamment de ce modele 
de modulation numerique lineaire : 

_ la modulation par codage a saut de phase (PSK pour "phase 
shift keying"), 

_ la modulation par codage a saut d' amplitude (ASK pour 
15 "amplitude shift keying"), et 

_ la modulation d'amplitude en quadrature (QAM pour -Quadra- 

ture Amplitude Modulation" ) . 

Avant d-en aborder la description detaillee, on enoncera 
20 d-abord les traitements proposes d'une maniere generale. 

A u niveau du principe, on estime d'abord les differentes 
inoonnues du modele. en supposant que le signal intercepts 
correspond bien a une modulation numerique lineaire; ensure, 
25 a partir de ces estimees. on validera ou non cette hypothese. 

Pour y arriver. la Demanderesse s'est interessee d'une part 
aux techniques de calcul matriciel sur des nombres complexes 
(algebre lineaire), d- autre part a des aspects statxstiques 

30 du traitement de signal. 

Les symboles contenus dans le signal de depart sent tir^s 
d *un -alphabet" de symboles a priori inconnu, avec une 
modulation elementaire (la modulation appliquee a chaque 
35 symbole) qui est egalement inconnue. 

Cependant, la Demanderesse a tout d'abord observe que le 
signal de depart etant de duree finie, le nombre de symboles 
successifs qu'il contient est egalement fini. L a Demanderesse 
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a alors pens<§ qu' apres avoir suppose que ce signal de depart 
verifie le modele de modulation numerique lineaire, il doit 
exister un moyen - au moins - de retrouver la suite des 
symboles, en s ' interessant au nombre de degres de liberte (ou 
parametres) qui sont necessaires pour obtenir une representa- 
tion du signal de depart (apres demodulation approximative 
dans 1 ' exemple d€crit ) . 

La Demanderesse a recherche un tel moyen. Ce faisant, elle 
s'est interessee tout d'abord a la maniere de representer le 
signal de depart (de preference apres une premiere demodula- 
tion). Comme on le verra ci-apres, elle a trouve des "bases 
de fonctions" qui permettent, par "projection", de represen- 
ter la plupart des signaux que l'on peut rencontrer en 
pratique, et ceci avec une perte d' informations minime, 
n^gligeable en l'espece. 

A partir de la, elle a developpe une technique, non limitati- 
ve, reposant sur des proprietes interessantes d'une solution 
particuliere du modele de modulation numerique lineaire; pour 
la determination des symboles, la Demanderesse s'appuie 
egalement sur le fait que ceux-ci sont decorreles ou faible- 
ment corr616s entre eux. 

Ladite solution particuliere, sur laquelle on reviendra ci- 
apres, est unique, et permet d'engendrer, par des combinai- 
sons lineaires, toutes les solutions possibles, comme on le 
verra egalement. II s'agit de la solution h^t) dont la 
longueur du support est la plus petite. 

On sait que, mathematiquement , le "support" d'une fonction 
f(x) peut etre d«§fini comme 1' ensemble des valeurs de la 
variable x pour lesquelles la fonction f (x) est differente de 
zero. On traite ici de fonctions du temps representant un 
signal de duree finie. Ainsi, le signal z(t) a un support 
borne; de meme, la fonction de modulation Slementaire modulee 
g m (t) a un support borne, que l'on peut appeler "longueur 
temporelle" de la modulation elementaire. De son cote, le 
support de la f onction- h^t) , c'est a dire la longueur de 



1-intervalle de temps sur lequel h^t) est non nulle. est au 
pl us egal a la longueur temporelle de la modulation eiemen- 
Lire on le comprend, du fait que la modulation elemental 
doit pouvoir se decomposer selon la fonction h„(t). 

on verra encore qu-apres application de cette solution h^t, 
d e support minimal, les symboles c m (k, assoc.es a cette 
solution particuliere s'expriment d'une maniere relat.vement 
simple en fonction des vrais symboles a«M du r.s.du de 

porteuse Af 0 - f. - f. * d * flltr "^ * 

permettront d'acceder a ces inconnues. 

II serait concevable de travailler sur la base d-une solution 
autre que celle de "support minimal", mais Impression du 
residu serait alors plus compliquee. 

L es dif f erents stades du traitement decrit en detail ci-apres 
sont (figure 3): 

_ passage en signal complexe z(t) (100), 

_ estimation du rythme 1/T (200), f 

. estimation -approve" «. de la frequence po rteus f (300 . 

- demodulation du signal complexe z(t) par cette tr q 
porteuse approchee f a (400), 

- estimation de la fonction h_<t> ( 500, 50! , 510, 520) , 

. estimation des symboles modules c»(k) assocres a h^t) 

l 61 esUmation du residu de filtrage a, et des symboles 
a partir des symboles c m (k) (650), 

_ estimation du residu de porteuse Af 0 < 700 >' DOrteuse 
0 - demodulation des symboles a.(k) par le residu de porteuse 
(600) afin qu'il ne reste que des symboles a(k), 
1 esllmLion des etats ou valeurs possibles e, des stoles 
a(k) a l'aide de l'alphabet (900), 

. estimation du train de symboles ef f ectivement em.s e(k) en 
55 fonction des etats possibles (950), n ^ ion 
- determination de 1'adequation du modele d'une -du^on 
numerique lineaire (990), utilisant les parametres estxmes 
precedemment (validation ou non de 1'hypothese) . 



Pour faciliter la comprehension, on decrira 1' ensemble des 
traitements, depuis la reception du signal eiectromagnetique. 
Mais il doit etre clair, comme deja indique, que la mise en 
oeuvre de 1' invention peut se faire a partir d'un signal 
enregistre, ayant deja subi une partie des traitements, en 
particulier le passage au signal complexe, auquel cas il est 
deja disponible en memoire 100. 

La figure 1 illustre le principe du passage en signal 
complexe (100). Le signal reel recu, tel qu'issu d'une 
antenne 10 se definit par 1' expression El. Pour les traite- 
ments, on lui fait subir un filtrage complexe 17, dont la 
fonction de transfert est adaptee au spectre du signal comme 
on le verra ci-apres. Ce filtrage fournit a un calculateur 20 
un signal complexe z(t) possedant une composante "en phase" 
I et une composante "en quadrature" Q. L'homme du metier sait 
que ce signal complexe peut etre represente comme les nombres 
complexes . 

Un filtre complexe centre sur une frequence f 1 peut etre 
defini par sa reponse impulsionnelle, qui correspond a la 
relation E2 : on part de la r6ponse impulsionnelle p(t) d'un 
filtre passe-bas reel classique (Fig. 4A) , que l'on multiplie 
par 1' expression E3, d'ou le spectre decaie de la figure 4B. 
Le mot "reel" signifie ici, comme dans la suite, "a valeurs 
r<5elles" , par opposition aux valeurs complexes. Ainsi, en 
reference a la figure 1, si le spectre du signal d' entree est 
celui de la figure 4C, un filtre ayant la fonction de trans- 
fert de la figure 4B donne un signal complexe ayant le 
spectre de la figure 4D. 

Dans un mode de realisation (figure 1A) , on prevoit, apres 
l'amplificateur de reception 11, un ou des Stages 13 a fre- 
quence intermediate (FI), suivis d'un Stage 15 de conversion 
en bande de base, puis du filtrage complexe 17. La conversion 
analogique numerique se fait generalement soit a 1' entree de 
l'etage 15, soit a 1' inter ieur de celui-ci, soit a sa sortie 
(avec une cadence d'echantillonnage plus basse dans ce 
dernier cas ) . 



D'autres techniques sont envisageables pour obtenir le signal 
complexe, en particulier une transformee de Hilbert (non 
representee), ou bien une demodulation complexe IQ en une 
composante en phase I et une composante en quadrature Q, 
comme illustre en figure IB s apres un etage 14, oil le signal 
se trouve dans une bande B autour d'une frequence centrale f 0 , 
on prevoit un Stage 18, avec deux voies de melangeurs 181 et 
182, recevant directement ou en quadrature (Stage 180) un 
signal cos 2rcf 0 t ; les mSlangeurs sont suivis de filtres 183 
et 185, qui donnent les sorties I et Q, respect ivement . La 
conversion analogique numerique se fait generalement avant ou 
apres 1' Stage 18. 

Le stade 200 consiste en une estimation du rythme 1/T de la 
modulation inconnue qui est prSsumee exister dans le signal 
recu. On decrira maintenant l'une des techniques utilisables 
a cet effet. 

A partir du signal complexe z(t), on forme en 201 (figure 6) 
0 le signal x(t) de l'expression E4, ou le decalage temporel r 
est petit devant la valeur T Sgale a 1' inverse du rythme que 
l'on cherche. On sait que, pour conjuguer z(t - r) dans E4, 
il suffit de changer le signe de sa composante imaginaire 
complexe. Le decalage temporel r peut etre choisi egal au pas 
5 d'Schantillonnage . 

La Demanderesse a relevS que le rythme recherchS correspond, 
dans le spectre de ce signal x(t), a une raie de frequence 
(1/T). Pour la suite des traitements, il est utile d'estxmer 
0 la position de cette raie, aussi precisement que possible. 

Pour cela, on calcule en 202 le spectre X(f) du signal x(t). 
Avantageusement, on utilise pour cela une transformee de 
Fourier rapide (FFT) entrelacee (technique dite du "zero 
15 padding", bien connue de l'homme de l'art). 

On choisit comme parametre de la FFT un entrelacement d'au 
moins 8 ( huit ) . 
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La sortie de cette FFT (Figure 5A) est ensuite traitee en 203, 
par exemple conf ornament h la relation E5 : on norme chaque 
valeur de X(f) obtenue & une frequence f par la moyenne des 
valeurs X(f) obtenues pour des frequences situees dans une 
5 fenetre de largeur 2.Af ± , centree autour de f . Ce traitement 
const itue une sorte de filtrage non lin6aire de la sortie de 
la FFT qui vise Si faire ressortir les raies contenues dans ce 
spectre (Figure 5B). La largeur 2Af ± fixe le degre de filtrage 
obtenu. Ici, elle a ete fix6e par simulation & 20 fois 
10 1' inverse de la duree totale du signal analyst. 

Ensuite, en 204, on seuille ce signal de fagon h seiectionner 
les raies qu'il contient, et, en 205 , les raies aux frequences 
negatives sont ignorees (Fig. 5C). 

15 

A cause de 1 ' entrelacement de la FFT, il reste plusieurs sous- 
raies au voisinage de f = 1/T (Fig. 5D) . La valeur 1/T 
recherchee est alors estimee a l'etape 206 par une interpola- 
tion qui peut consister h ajuster une courbe, du type 
20 polynomiale, par exemple, qui passe par les sommets de ces 
differentes sous-raies et en prenant le maximum de cette 
courbe . 

En paralieie sur 1' estimation du rythme de la modulation, on 
25 determine une estimee approchee f a de la frequence porteuse 
fQ (stade 300). L'estimee approchee est utile, dans la mesure 
ou elle servira & effectuer une premiere demodulation, qui 
permet de travailler ensuite & basse frequence : la base de 
fonctions sera definie en basse frequence, et les calculs 
30 matriciels seront effectu6s sur des echantillons egalement 
disponibles h basse frequence. 

L'estimee approchee f & peut etre obtenue de differentes 
maniferes. La Demanderesse prend de preference pour f a la 
35 valeur moyenne de la frequence instantanee du signal. Cette 
valeur moyenne f a est par exemple calcuiee comme suit: 
- comme en 201, on prend d'abord le signal x(t) de 1' expres- 
sion E4, puis 
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- en 311, on calcule la valeur moyenne de la phase de x(t), 
sur ensemble de la duree du signal, et on la divise par 

5 Cette technique repose sur le fait que, si r est petit devant 
T l'ecart entre la phase du signal z(t) a 1' instant t et la 
phase a 1' instant (t-t) ne depend que de r et f 0 . Comme 
precedemment, r peut etre choisi egal au pas d'echantillon- 
nage (cas de la figure 6), ou bien calcule en prenant une 

10 fraction de T petite „ 

La figure 7A montre un exemple, ou x(t) est definie par un 
nuage de points, correspondant a une valeur moyenne f a - 
1940 Hz. 

Si 1' on a estime f a , on precede (stade 400) a une demodulation 
du signal complexe z(t) par cette estimee approchee f a de la 
porteuse f 0 , de maniere connue. Numeriquement , la demodulation 
consiste en une multiplication complexe du signal z(t) par 
une expression semblable a E3, ou f, est remplacee par -f a . 

Apres cette demodulation, le spectre initial (Fig. 7B) du 
signal complexe vient se centrer approximativement sur la 
frequence zero (Fig. 7C). Le decentrage residuel du spectre 
correspond a un petit residu de porteuse Af 0 Sgal a f 0 - fa- 
un exemple est donne en figure 7D, ou le trait mixte re P r6- 
sente le spectre du signal avant demodulation, tandis que le 
trait continu represente le spectre du signal apres demodula- 
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On notera a nouveau z(t) le signal obtenu apres demodulation 
par f . Conf ormement au modele de modulation numerique lmeai- 
re ce signal z(t) satisfait desormais la formule E10. La 
Demanderesse a observe que ce signal z(t) peut aussi s'ecrxre 
selon 1' expression E14, comme une somme de produits dont 
chacun fait intervenir: 

- une fonction de modulation Element aire modulee g m (t), 
d6finie par la relation E12, et 



- une suite de symboles modules a^k), d6finie par la relation 
E13. 

Dans ces expressions et par la suite, l'indice m signifie 
"module" , et k va de 1 a kx, qui est le nombre total de 
symboles intervenant dans le signal analyst. 

La demonstration de ceci passe par la r6-ecriture de la 
formule E10 de la facon exprim^e dans la formule Ell. 

On notera que l'expression E14 est "sans porteuse apparente" , 
c'est a dire que le r<§sidu de porteuse Af 0 n'y apparait plus 
explicitement. En effet, il est contenu dans la fonction de 
modulation eiementaire "modulee" g m (t) et dans les symboles 
modules a^k) . 

Pour acc£der aux symboles, il y aura lieu d'estimer cette 
fonction de modulation eiementaire g m (t). Sont 6galement 
inconnus le r£sidu de porteuse Af 0 , 1' ensemble des symboles 
utilises, ou "alphabet", et la suite ou train de symboles 
precisement contenus dans le signal de depart. 

Dans le mode de realisation decrit, le stade 500 (501-520) 
consiste en 1' estimation de la fonction de modulation 
<51<§mentaire modulee, c'est a dire de la fonction g m (t). On en 
deduira les symboles modules a^k) au stade 600 (610-650) 
(voir figure 3). 

Se presente d'abord la difficulty suivante : travailler avec 
un signal defini pour chaque valeur de t est une tache 
difficile, compte-tenu du nombre potentiellement tres eieve 
de valeurs a considerer. 

La Demanderesse propose de proceder a une approximation, en 
projetant le signal z(t) sur une « base » de fonctions, 
choisie a l'avance. Cette base (au sens mathematique du mot) 
est en fait une famille libre, de preference orthonorm6e, de 
fonctions ( ip x ( t ) , <p 2 ( t ) , . . . , tp h ( t ) ) . 
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On choisit une base de fonctions dp lf *2r • • -*b> adaptee au 
rythme 1/T tel qu'estime precedemment , et construite par 
exemple de la facon indiquee dans les expressions E25 et E26. 
Dans ces expressions, rect< ) designe la fonction rectangu- 
laire sur 1'intervalle de temps indique en dessous. L'entxer 
q est tel que le produit q.T soit sensiblement egal a la 
longueur du signal analyse, b correspond au nombre de 
fonctions de base necessaire pour couvrir toute la duree du 
signal analyst- 

D'autres bases de fonctions envisageables seront indiquees 
plus loin. Toutefois, la Demanderesse considere actuellement 
aue la base de fonctions des expressions E25 et E26 demeure 
la plus simple, et donne des approximations largement 
suffisantes pour les signaux rencontres en pratique. 

Quel que soit le mode de realisation choisi, il existe dans 
le dispositif une "representation" de la base de fonctxons. 
C-est ce qui est rappele par le bloc 501 de la figure 3. En 
1'espece, cette representation est parametrable par T. 

Dans le mode de realisation decrit en detail ci-apres on 
utilise la base definie par E25 et E26. On a ainsi deux 
fonctions de base par intervalle de T secondes , on obtiendra 
ainsi deux coefficients de projection toutes les T secondes 
( on dit aussi "deux points par symbole", ttant donne que 1'on 
recoit un symbole toutes les T secondes). On verra plus loxn 
comment 1' invention peut s'etendre a des valeurs entieres 
super ieures h 2 (q^2) * 

L ' approximation commence alors par le calcul de la projection 
du signal z(t) dans L'eapace (egalement au sens mathematxque 
du mot) engendre par les fonctions ^(t). Ce signal P ro,ete 
ne sera alors plus defini que par des coefficients de projec- 
tion, qui eux, seront en nombre fini, inferieur au nombre 
d'echantillons temporels du signal z(t). 

Le signal z(t) est done tout d'abord projete (510, fig. 3) 
sur la base de fonctions definie par les relations E25 et 
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E26. Cette projection est notee z p (t), ou l'indice p signifie 
« projet6 » . 

La Demanderesse a observe qu ' approximer le signal z(t) par sa 
projection dans la base des ^(t) revient a approximer les 
fonctions g m (t-kT) (k variant) intervenant dans le modele du 
signal par leurs projections dans cette meme "base". On 
choisit la base de telle sorte que toutes les fonctions de 
base se d«§duisent par translation d'un nombre entier de T 
d'un petit nombre de fonctions <p 1 (t) , v>2< t > ' • • • ( deux f onc ~ 
tions pour une base a deux points par symboles, q fonctions 
dans le cas general). La Demanderesse a alors observe que les 
projections des fonctions g m (t-kT) se d<§duisent les unes des 
autres par translation de T: elles sont done dSfinies par les 
memes coefficients de projection a un decalage pres dans le 
temps. Tout se passe comme si on avait approxinul la fonction 
inconnue g m (t) par sa projection dans la base ^(t), qui 
s'ecrit avec un petit nombre de coefficients inconnus, qu'il 
s'agira de trouver. 

Pour plus de simplicity, on choisit les fonctions ^(t) , 
p 2 <t),... generant la base par translation de kT (k e N) , de 
telle sorte qu' elles se d<§duisent, elles aussi, par transla- 
tion de T/2 (plus generalement de T/q, q e N) d'une meme 
fonction ^ ± (t). On obtient ainsi les bases E25 et E26 (deux 
points par symbole) et E34 et E35 (q points par symbole) , ou 
les fonctions ^(t) se d«§duisent finalement les unes des 
autres par translation respectivement de kT/2 et kT/4 (k e 
Z) . 

La Demanderesse a observe en outre que les fonctions g m (t) 
rencontrees en pratique sont suffisamment lisses pour etre 
approximees avec peu d'erreurs dans ces bases de fonctions. 
Un exemple schematique d'une telle approximation est illustr6 
sur la figure 8C. 

La Demanderesse a aussi observe que, pour projeter le signal 
z(t) dans la base de fonctions dSfinie par les relations E25 
et E2 6, il est possible de filtrer le signal par le filtre de 

0( 



reponse impulsions lie ^(-t) (expression E25), puis d'echan- 
tillonner la sortie de ce filtre toutes les T/2 secondes. II 
apparait done que la realisation de cette projection est 
particulierement simple. 

Comme le montre la figure 10A, le signal est centre sur la 
frequence zero et possede, dans les cas pratiques , une bande 
spectrale de largeur 1/T. Le filtre de projection est aussi 
centra sur la frequence zero ; par centre, il possede une 
bande de 2/T (a 4 dB) . L'homme du metier en d<§duira, comme 
cela a d6ja 6t<§ dit pour la fonction g m (t), que 1' approxima- 
tion faite par cette projection a pour effet que la partie du 
signal perdue par projection est negligeable. 

incidemment, on remarquera des synergies interessantes de ce 
mode de realisation: 1' estimation prealable du rythme permet 
que la projection du signal sur la base de fonctions soit 
bien synchronised sur le rythme de modulation ; 1' estimation 
de la moyenne de la frequence instantanee permet que cette 
projection se fasse sur une base de fonctions de basse fr6- 
quence * 

Les travaux de la Demanderesse ont montr<§ qu'avec deux 
<§chantillons par symbole (au moins), e'est-a-dire avec au 
moins deux fonctions de projection toutes les T secondes, il 
devient possible de trouver des solutions du modele (de 
modulation num^rique lineaire) qui correspondent au signal de 
depart, a l'aide de techniques de calcul matriciel (bloc 520, 
fig. 3 ) . 

Cependant, lesdits travaux ont montr# aussi que, pour un 
signal de depart donne, la relation E14 admet une multi- 
plicity de couples de solutions possibles. 

Parmi toutes ces solutions possibles, la Demanderesse s'est 
attached a rechercher celle qui verifie la propriete suivan- 
tes il s'agit de la solution donnant des symboles modules 
a m (k) qui, apres demodulation, prendront un nombre fini, et 
le plus petit possible, de valeurs. 

d 



AprSs avoir envisage d'imposer cette contrainte directement , 
la Demanderesse pr6ffere actuellement proc«§der en deux temps, 
comme on le verra. 

Au moins lorsque z(t) est un signal issu d'une modulation 
num^rique lin^aire, la Demanderesse a montr6 ce qui suit : 

- parmi les nombreuses solutions possibles de la mod^lisa- 
tion, il n'en existe qu'une seule (c m (k) ,h m (t ) ) v^rifiant la 
relation E20 f ou h m est, parmi les fonctions de modulation 
solutions, celle dont la longueur du support est la plus 
petite ; 

- en outre, cette solution particulidre (c m (k) ^(t) ) engendre 
1' ensemble des solutions possibles de la fagon suivante: si 
(a m (k) ,g m (t ) ) est une autre solution, alors il existe n 
coefficients (a x , . . .a n ) permettant de verifier 1'expression 
E21 donnant g m (t), tandis que les symboles a^k) sont d^finis 
r^cursivement, pour tous les k, par 1'expression E23. 

Elle en a tir§ que la fonction de modulation 616mentaire 
modul^e g m (t) peut s'<§crire selon 1'expression E21, ou h^t) 
est 1' unique fonction solution de support le plus petit. 

S'appuyant sur les remarques qui precedent, la Demanderesse 
pr<§conise de decomposer la recherche de la fonction g m en deux 
parties : 

dans un premier temps, on recherche la solution 
(c m (k) ,h m (t) ) avec la contrainte que h m (t) ait le support le 
plus petit possible, 

- ensuite, & partir des symboles c m (k) correspondant h cette 
solution h^ et se mod61isant selon 1'expression E2 3, on 
estime la deuxieme partie de la fonction g m , h savoir les 
coefficients (a ± ) , qu'on nomme globalement "r6sidu de fil- 
trage" . 

On mod61isera par la suite le signal par le couple de 
relations E22 et E23, ou h^t) est la fonction solution de 
support le plus petit, et les (a ± ) sont le rSsidu de filtrage. 
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Comme cela a ete dit precedemment , approximer le sxgnal z(t) 
par sa projection dans 1'espace des fonctions ^(t) revxent 
a approximer la fonction g m (t) par sa projection dans ce m£me 
espace. II en est de meme pour la solution particuliere h^t), 
que 1'on va remplacer par sa projection dans la base des 

fonctions «Mt) ,*> 2 (t) *>(*>>• Cette Projection h pm (t) ne 

sera alors plus definie que par un petit nombre de coeffx- 
cients de projection . 

) Maintenant, la projection z p (t) de z(t) ne se modelise plus 
qu'a partir de la projection h pm (t) de h^t). La (petxte) 
partie du signal qui est orthogonale a la « base » de 
projection est perdue. La base de fonctions E25, E26 a ete 
choisie de telle sorte que la partie du signal perdue soxt 

5 suffisamment faible pour les signaux rencontres en pratxque. 

On donnera maintenant, en reference a la figure 9, un exemple 
du deroulement des traitements visant a estimer h^t). 

0 Tout d'abord, comme deja indique, le signal z(t). est projete 
sur la base de fonctions <p^) . En 1'espece, cela consxste, 
en 510, a filtrer le signal z(t) par le filtre de reponse 
impulsionnelle ^(-t) (expression E25 ) , puis, en 523, a echan- 
tillonner la sortie de ce filtre toutes les T/2 secondes, 

25 jusqu'Si la fin du signal z(t). 

On estime ensuite la fonction h pm (t) de support le plus petit, 
qui soit solution de 1' equation E30. 

La projection z p (t) du signal z(t) peut s'ecrire a partir de 
ses coefficients de projection z(k) selon la relatxon E36 
(avec p pour projet^)- 

De son cote, la projection h pm (t) de la fonction de support 
le plus petit h.(t) s'ecrit selon la relation E37 ou L est 
la longueur de h pm (t), comptee en nombre de perxodes ^ T . 
Autrement dit, la longueur de h pm (t) est superieure a (L-l).T 
et infer ieure ou 6gale a L.T. 
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Compte-tenu de la cadence de deux dchantillons (deux coeffi- 
cients de projection) par symbole toutes les T secondes, deux 
gchantillons successifs de rang 2k + 1 et 2k + 2 de la 
projection z p (t) correspondent au couple de valeurs z(2k+l) 
et z(2k+2), lequel couple se mod61ise selon les relations E38. 

La premifere 6tape du traitement consiste a estimer la 
longueur L de h pm (t). 

A cet 6gard, la Demanderesse a remarqu6 que les longueurs des 
fonctions de modulation «§l<§mentaire g m (t) rencontr6es en 
pratique sont g6n6ralement inf^rieures & 10 fois T. 

On choisit alors une borne sup^rieure param^trable L lf telle 
que L x = 10.T. A chaque rang k, l'<§tape 525 fait correspondre 
un vecteur ou matrice unicolonne Z L1 (k) conforme & 1' expres- 
sion E41. Les vecteurs Z L1 (k) correspondent Sl des tranches de 
signal z(t) prises entre des instants kT et (k + L1)T. Ce 
vecteur Z L1 (k) contient 2.L1 composantes successives en r(k), 
allant de (2.k+l) & (2.k+2.Ll). Les vecteurs Z L1 (k) se 
recouvrent en grande partie les uns les autres: Z L1 (k+l) se 
deduit de Z L1 (k) en supprimant en haut z(2k+l) et z(2k+2), en 
faisant remonter les autres elements de deux cases, et en 
ajoutant en bas z(2k+2Ll+l) et z(2k+2Ll+2). La construction 
de ces vecteurs Z L1 (k) & partir des donnSes de base se fait 
done avec un entrelacement . 

En 527 f on forme la matrice carr6e M Llf de dimension 2. LI, 
d<§finie par la relation E42. II s'agit d'une sorte de matrice 
de covariance des grandeurs z(k); cependant, k la difference 
de la covariance habituelle, les vecteurs constitutifs sont 
entrelaces ; e'est pourquoi la Demanderesse propose de nommer 
M L1 la "matrice de covariance entrelac^e" . 

On calcule en 52 9 les valeurs propres de M L1 , de mani&re 
connue. L' analyse de celles-ci permet d'acc^der h la longueur 
L de la fonction h pm (t) f car la Demanderesse a 6tabli que la 
matrice M L1 poss^de L x + L - 1 valeurs propres positives non 
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null... t L t - L + 1 (infSrieur a L, + L - 1) valeurs propres 
null.., en fait presque nulles a cause du bruxt. 

Par exemple, on P eut ranger ces valeurs propres par ordre 
decroissant de grandeur. La repartition des valeurs propres 
est approximativement donnee par la figure 12, ou 1'ordonnee 
est la grandeur des valeurs propres, lesquelles sont rangees 
par ordre decroissant de grandeur en abscisse. On estxme 
1'indice ou rang de la valeur propre (apres tri) ^f^' 
pond au point de cassure de la ligne de la figure 12, c est 
I dire au point ou 1'on descend brusquement a des valeurs 
propres presque nulles, correspondant au bruit Ce point de 
cassure est estime classiquement en utilisant le 
"MDL" (pour Minimum Description Lenght). On pourra se referer 
a 1'article "Detection of Signals by information Theoretic 
criteria", IEEE Transactions on Acoustics, Speech and Signal 
Processing, Avril 1985, P 387. L ■ estimation de ce point ^ de 
cassure peut etre interprets comme 1' estimation du rang de la 
matrice 

La longueur L de la fonction h^t) est donnee par le rang 
apres tri de la premiere valeur propre nulle. 

on forme maintenant en 535 des vecteurs z L <k, definis par les 
2.L composantes consecutives en .(k) de 1-expressxon E44. En 
537, on forme. I partir de ces veoteurs M*., la 
(dimension 2L x 2L> de 1'expression E43 (sxmxlarre a 1 ex 
pression E42, mais avec L au lieu de Lj). 

cette matrioe est une autre "matrice de covariance entrela- 
cee», mais avec un entrelacement restreint a I. au lieu de LI. 
Les veoteurs z L (k) correspondent a des tranches du signal 
« (t) prises entre les instants kT et (k + L)T. 

Le traitement est alors le suivant. En 539 on determine la 
plus faible valeur propre de M L en calculant, par exemple 
comme precedemment , toutes les valeurs propres (qui. sont 
positives a cause de la symetrie hermitienne de M L ) et en les 
rangeant par ordre decroissant. 
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En 541 , on calcule le vecteur propre associe h cette valeur 
propre la plus faible. La Demanderesse a en effet 6tabli: 

- que la matrice M L a dans ce cas une seule valeur propre A 
presque nulle, 

- que la grandeur de cette valeur propre permet d'estimer la 
puissance du bruit o 2 = A / L, et 

- que le vecteur propre associ§ & cette valeur propre (ex- 
pression E45) ne depend que de la fonction h pm (t), solution 
ayant le support le plus petit. 

En 543, les coefficients de projection de h pm (t) sont done 
obtenus (au signe prfes) en prenant dans l'ordre les conjugu6s 
complexes des composantes de ce vecteur propre, et en 
changeant une fois sur deux le signe de ces composantes (on 
peut travailler au signe pres). 

On notera au passage que les 6tapes 535 a 539 sont les memes 
que 525 h 529, en remplagant LI par L. Elles peuvent done 
etre mises en oeuvre par les m§mes moyens de calcul. 

Ce traitement peut etre interpret^ comme suit : 

- les vecteurs Z L1 (k) correspondent a un d^coupage du signal 
z (t) en tranches de longueur L 1 .T, qui se recouvrent; 

- ces vecteurs sont d^finis par les 2.L X coefficients de 
projection (z(2k + 1 ),..., z( 2k + 2L1)), qui sont a priori 
susceptibles de prendre n'importe quelle valeur dans l'espace 
vector iel de dimension 2.L X correspondant ; 

- cependant, sur chaque tranche ce signal z p (t) devient le 
signal u(t) conforme & la relation E4 7, dont les 616ments 
z x (t) , . . . r^+jsi-iit) sont d6finis par les relations E48; 

- par consequent, sur chacune de ces tranches, le signal z p (t) 
est tou jours une combinaison lin^aire des m§mes (L x + L - 1) 
signaux z x (t ),..., z L +Li-i( t ) ' combinaison lin6aire dont les 
coefficients sont les symboles c m (k) intervenant dans la 
tranche consid6r6e et variant done d'une tranche h 1 'autre ; 

- bien que d^fini au depart h partir de ses 2.L X coefficients 
de projection (z(2k + 1 ),..., z( 2k + 2L1)), chaque tranche du 
signal z p (t) ne prend done en fait ses valeurs que dans le 
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meme sous-espace de fonctions engendre par les L x + L - 1 
signaux definis ci-dessus. 

La Demanderesse a montre aussi que la condition « h pm est la 
solution ayant le support le plus petit » est equivalente a 
la condition « les L 1+ L-l signaux precedents sont lineairement 
independants » . 

Sur chaque tranche de longueur L X T, le signal z p (t) prend done 
ses valeurs dans un espace de dimension L x +L-1. 

La matrice M L1 definie precedemment permet, par diagonali- 
sation, d'estimer le sous-espace contenant ces tranches de 
signal, a partir des vecteurs propres associes aux valeurs 
propres non nulles. Pour simplifier le traitement, on peut, 
comme decrit, n'estimer que le rang L 1+ L-l de cette matrice 
M L1 (nombre de valeurs propres qui ne sont pas tres faibles 
ou nulles). Apres en avoir deduit L, il est alors possible de 
redecouper le signal en tranches de longueur L.T. 

Le meme raisonnement que precedemment permet alors de voir 
que ces tranches de signal engendrent, dans I'espace de 
fonctions de dimension 2.L, un sous-espace de dimension 2.L- 
1, car: 

- on sait maintenant que chacune des tranches de longueur L.T 
du signal Zp (t) est susceptible a priori de prendre ses 
valeurs dans un espace de dimension 2.L mais qu'en fait 

- sur chacune de ces tranches, le signal z p (t) est toujours 
une combinaison lineaire des memes (2.L-1) signaux 
s^tO.-.B^iCt) de 1'equation E48 ou l'on remplace Ll par L. 

En bref , il existe, dans I'espace de fonctions, une direction 
qui est orthogonale a chacune des tranches de signal s e'est 
celle engendree par le vecteur propre V de 1' expression E45, 
35 associe a la valeur propre nulle de la matrice M L . 

On a maintenant estime h pm (t). On peut alors passer au stade 
600, qui est l'estimation du train de symboles modules c m (k), 
associes a cette fonction h pm (t) estimee. 

0/ 
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Ceci se fait en filtrant le signal z p (t) par le filtre inverse 
610 de h pm (t). On indiquera main-tenant une maniere, non 
limitative , de calculer le filtre inverse, en passant par la 
matrice pseudo-inverse de la matrice de filtrage associ^e k 

On part de 1' expression de h pm (t) donn^e en E58 qui est le 
r£sultat du traitement pr6c6dent, obtenue par son extension 
pour une projection & plus de deux points toutes les T 
secondes. Pour le cas particulier de deux points par symbole, 
traits jusqu'& present, il suffit de poser q = 2. 

On forme ensuite la matrice H, donn6e par 1' expression E63. 
II s'agit d'une matrice rectangulaire de dimension qL. (LI + 
L-l) que l'on peut appeler "matrice de filtrage". Elle est 
form«§e de la fagon~ suivante. Chaque colonne possfede q.L 
composantes. La premiere colonne est obtenue en prenant les 
q derniferes composantes de h pm (t), et en compliant avec des 
z«§ros. Pour le reste, chaque colonne se d6duit de la pr6c6- 
dente par un d6calage de q cases vers le bas, en compliant 
en haut par les composantes de h pm (t) dans l'ordre, puis par 
des z<§ros, lorsque ces composantes sont £puis6es. On continue 
ainsi jusqu'a la derniere colonne ne contenant en haut que 
des z«§ros puis en bas les q premieres de h pm (t). On obtient 
ainsi une matrice avec LI + L -1 colonnes. 

On forme alors le produit matriciel de 1' expression E64, ou 
l'homme du metier reconnaitra la matrice pseudo-inverse de H, 
de dimension (2.L-l).qL. Les coefficients du filtre inverse 
sont estim«§s par la ligne numSro L de la matrice donn^e en 
E64. 

On applique alors en 610 le signal z p (t) au filtre inverse 
ainsi obtenu. La sortie correspond aux symboles c m (k) de 
1' expression E43, dans laquelle a ± est le r6sidu de filtrage 
et les a^k) correspondent aux symboles recherch^s, modules 
par le r^sidu de porteuse Af 0 , toutes ces grandeurs 6tant 
encore inconnues h. ce stade de 1' algorithms. 



Pour estimer ces inconnues, comme on 1'a vu, la Demanderesse 
a envisage d'imposer la condition suivante : les symboles 
a-(k), associes au residu de filtrage (a L ) . doivent correspon- 
dre a des symboles modules, qui, apres demodulation, pren- 
dront un nombre fini, et le plus petit possible, de valeurs. 
La condition que la Demanderesse prSfere imposer aujourd'hux 
est • "le train de symboles (a.(k)) est, parmi les solutions 
possibles donnant des symboles decorrel<5s, celui pour lequel 
1' alphabet (ensemble des valeurs possibles) est de taxllo la 
plus petite" . La condition de decollation des symboles est 
donnee par la relation E70 et porte sur la covariance des 
symboles. La deuxieme condition consists en fart a chdsxr 
parmi les solutions permettant de surmonter la premiere 
contrainte, le train de symboles dont la variance du module 
des symboles est la plus faible (equation E80). 

Ces deux contraintes, plus simples a mettre en oeuvre, 
aboutissent sensiblement au meme resultat pour les sxgnaux 
rencontres en pratique. (Les vrais symboles ends sont 
g6n6ralement d6corr61es.) 

L'obtention des elements individuels (a ± ) du residu de 
filtrage s'effectue en 650. 

A 1'etape 651 (figure 15), la Demanderesse part d'un polyn6me 
Q(Z) construit a partir des mesures conformement a 1' expres- 
sion E73, ou Z est la variable. Ensuites 

- en 653, on calcule les N racines complexes fi L de module 
infer ieur a 1 de ce polynSme Q(Z), 

- en 655, on construit les 2* possibility differentes de N- 
uplets construites a partir des N racines p t . en remplacant, 
ou non, chaque racine fi ± par son inverse conjugue l/fi L . 
Chacun des N-uplets est alors interprets comme les racines 
d'un f iltre, 

- en 657, on obtient ainsi 2 N f litres possibles dont 1 un 
correspond au r<§sidu de filtrage a L que 1'on cherche a calcu- 
ler, 

- puis, en 658, pour chaque possibility de f iltre ainsi 
obtenue, on estime a partir des symboles c m (k) le tram de 
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symboles a^k) qui lui correspond (par filtrage inverse) 
conf ornament & la relation E23. 

II reste alors en 659 h choisir, parmi ces 2 N trains de 
symboles possibles, celui qui correspond h des symboles a m (k), 
qui, aprfes demodulation par le residu de porteuse Af 0 encore 
inconnu h ce stade, prendront le nombre, le plus petit 
possible, de valeurs dif f erentes . Une fagon d'operer consiste 
alors k determiner, parmi tous les trains de symboles 
possibles ainsi identifies, celui dont la variance moyenne du 
module des symboles (encore modules) est la plus petite. 
Cette variance est definie par la relation E80, dans laquelle 
nb designe le nombre de symboles. 

La base de ce traitement va maintenant etre expliquee. 

Au residu de filtrage (a ± ), on peut associer un polyn6me P 
correspondant & sa transformee en 2 dont la forme est decrite 
par 1' expression E71, ou N est un majorant du nombre n de 
coefficients a ± . En effet, ce nombre n n'est pas connu h 
l'avance, mais on sait, pour les signaux rencontres en prati- 
que, lui trouver un majorant N. On a dej& vu, en effet, que 
les formes d'impulsion eiementaires g(t) rencontres en 
pratique ont une longueur infer ieure h 10.T. D'apres 1' ex- 
pression E12, il en est done de meme pour g m (t) . En utilisant 
1' expression E21, on en deduit done que m est inferieur & 10 
( la longueur de h^t) est en nombre de T au moins egale & 1) * 
On prend en pratique N = 10 comme majorant, mais cette valeur 
est parametrable dans 1 ' algorithme . Prendre un majorant peut 
seulement avoir pour effet que les derniers coefficients a ± 
avec i e {n+l,...,N} soient nuls, si ce majorant N s'avere 
etre super ieur a n. 

On definit en plus le polynome P 2 (Z) de 1' expression E72, 
correspondant au module carre du polynome P(Z), et que l'on 
peut interpreter comme la transformee en Z du module carre de 
la fonction de transfert du residu de filtrage a ± . II est 
possible aussi d'ecrire les polyndmes P(Z) et P 2 (Z) sous forme 
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factorial, selon les expressions E75 et E76, oil les fi ± sont 
les racines du polynome P- 

La oemanderesse a montre, en utilisant 1'hypothese de decor- 
; relation des symboles, que le P olyn6me P 2 (Z) peut etre estxme 
par le polynome Q(Z) precis, a une constante multiplicative 
pres. Ce polynome Q(Z) peut etre, lui, interprets comme une 
approximation de la transformee en Z de la fonction d' auto- 
correlation des symboles. Bile a egalement observe que : 
D - si fi ± est une racine de P 2 alors 1//J ± * est aussi une racme 

^sl 2 ^ est une racine de P 2 alors fi L ou bien l/fi L * est une 
racine ^u polyndme P qu'on cherche a estimer. 

5 Dans une variante, il est possible d'etendre ce traitement au 
cas de symboles qui seraient legerement corrSles, c'est-a- 
dire eventuellement correles mais seulement avec leurs 
proches voisins, conformant a la relation E81. Dans ce cas, 
la Demanderesse a observe qu'une racine fi L de Q(Z) ne sera pas 

tO tou jours une racine de P 2 (Z), qui sert a definir les racxnes 
de P(Z) . 

Pour former les differentes possibility pour P, pour chaque 
racine fi L de Q et de module inferieur a 1 , on aura nontenant 

25 le choix entre prendre ou bien prendre l/fi L , ou bxen ne 

prendre ni fi L ni Si N est le degre du polyndme Q, on 

aura done 3»' f litres possibles pour choisir le resxdu de 
filtrage. Le choix de la bonne solution se fait alors comme 
precedemment, en minimisant la variance du module des 

30 symboles. 

On examinera maintenant un exemple de resultats de 1' ensemble 
des traitements d'estimation de g pm (t): 



- la figure 14A illustre la fonction h 



35 - la figure 14B illustre le residu de filtrage, ou se 
detachent des pics correspondant aux coefficients a L ; 
- la figure 14C illustre la fonction Slementaire de modula- 
tion modul<§e g pm (t) ; et 

d 



- la figure 17 illustre les symboles a^k) en sortie du 
traitement, qui correspondent aux symboles recherches, mais 
demeurant modules par le residu de porteuse Af Q . 



Ces symboles se modeiisent par 1' expression E13. II convient 
maintenant de trouver le residu de porteuse Af Q pour en 
deduire les vrais symboles a(k), entierement d§modules. Le 
traitement 700 consiste a tester tout ou partie d'une plage 
de frequences [-f lf fj, discretisee avec un pas suffisamment 
fin. Pour chaque valeur testee du residu de porteuse, on 
calcule la density de probability d(x,y) des symboles 
correspondants selon la formule E82, ou bien une autre 
grandeur de memes proprietes. On pourra se referer a 1' arti- 
cle "an estimation of a probability density function and 
mode", E.Parzen, Annals of Mathematical Statistics, 1962. Le 
critere utilise pour seiectionner la bonne valeur consiste a 
choisir la frequence Af Q telle que 1' expression E85 soit 
maximum. On peut assimiler la maximisation de ce critere a 
une minimisation de l'entropie (mesure du desordre) des 
symboles obtenues apres demodulation. II serait d'ailleurs 
possible de remplacer ce critere par d'autres fonctions de 
1 ' entropie . 

Quand la frequence testee est proche de la bonne valeur 
(Figure 16), la repartition dans le plan complexe des 
symboles a(k), qui s'en deduit par demodulation par la 
frequence testee, correspond a une constellation a nombre 
d'etats petit, d'ou une densite de probabilite " concent r6e " , 
mieux percue dans la representation tri-dimensionnelle de la 
figure 18. 

Quand la frequence testee est loin de la bonne valeur, la 
repartition des symboles a(k) correspondante est plus etaiee 
suivant des cercles centres sur l'origine (Figure 17, ou 
l'echelle deforme les cercles en ellipses). La figure 19 le 
montre en vue tri-dimensionnelle. 

Une fois la valeur de Af 0 ainsi estimee, une demodulation 800 
permet d'obtenir le train de symboles a(k) qui lui corres- 
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pond. Cette demodulation consiste en une multiplication 
complexe du signal par une expression semblable a E3, ou f x 
est remplac^e par -Af Q . 

5 II faut encore estimer les etats possibles des symboles 
(stade 900) . 

Ce stade tient compte du fait que les symboles obtenus apres 
demodulation correspondent aux vrais symboles, a nombre 
d' etats finis, mais corrompus par le bruit de puissance a , 
calculee precedemment , comme le montre la figure 16 qux 
correspond a un exemple de resultat obtenu. Les differentes 
valeurs des symboles a(k) trouv<5s sont d^ssinees dans le plan 
complexe. La figure 2D correspond a 1' alphabet des symboles 
a(k) reellement emis pour cet exemple. La rotation de phase 
entre les deux figures ainsi que les echelles differentes, 
proviennent du fait que, tout <§tant inconnu au depart, xl 
n'est possible de retrouver le train de symboles qu'a une 
constante complexe multiplicative pres. 
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L' object if de cette partie du traitement est d' estimer la 
position des etats possibles. Le deroulement du traitement 
est le suivant (figure 20): 

- on estime le premier <§tat possible ei comme la position en 
25 (x,y) du maximum de la density de probability d(x,y) des 

symboles bruitSs, 

- a partir de cette densite, on construit la fonction d 2 (x,y) 
correspondant a d(x,y), mais dont on force la valeur a zero 
pour tous les points se trouvant dans un rayon de 3o (par 

30 exemple) autour de e lf * 

_ on estime alors le deuxieme <§tat possible e 2 en prenant la 
position en (x,y) ou la fonction d 2 est maximale, 

- on continue ainsi le traitement jusqu'a ce que tous les 
etats possibles soient trait^s. 

Ceci peut se faire selon les Stapes 901 a 909 de la figure 
20. Parmi differentes possibility , le critere d'arret de ce 
traitement (condition "COND" en 907), permettant de savoir si 
tous les etats possibles ont M estim^s, peut etre celux de 
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1' expression E86, ou A est un seuil qu'on se fixe par avance 
(0 < A < 1) . 

Le principe utilise est le suivant. On suppose que les 
5 valeurs des symboles 6mis sont 6quiprobablement distributes 
entre les difftrentes valeurs possibles de 1' alphabet (on a 
& peu pr&s le meme nombre de symboles dans chacun des modes 
de la figure 18, par exemple). Les inttgrales de density de 
1' expression E86 sont proportionnelles aux nombres de 

10 symboles formant ces densitts . Supposons qu'on ait estim6 n - 
1 etats e lf ..., e n-1 . Le crit^re d'arret consiste h calculer 
une mesure du nombre moyen par 6tat de symboles correspondant 
aux 6tats e lf ..., e n-1 (membre de droite de E86) et h comparer 
cette quantite a la mesure du nombre de symboles restants. Si 

15 cette derniere quantity est trop petite, c'est qu'il ne reste 
pas assez de symboles pour former un etat suppltmentaire . Ces 
symboles sont supposes etre loins des 6tats e lf ..-, e n 
estim6s, & cause des pics de bruit. 

2 0 On peut maintenant proctder h 1' estimation de la suite des 
symboles, tel qu'elle a <§t6 ef f ectivement 6mise (stade 950). 

Cette partie du traitement consiste, & partir du train de 
symboles bruitts a(k) et des etats possibles des symboles e lr 
25 e 2 ,... calculus pr«§c<§demment , h estimer, pour chaque symbole 
bruits a(k), l'etat e(k) intentionnellement 6mis . 

A cet effet, pour chaque symbole a(k), en partant de k = 1 
(etape 951, figure 21) on recherche l'6tat possible e k le plus 

30 proche, par exemple (953) sur la base de la distance eucli- 
dienne (fonction DIST[ ] ) . On obtiendra ainsi, pour le rang k, 
un symbole not€ e(k) prenant une des valeurs e ± de 1' alphabet 
estimt pr6c6demment . Le minimum de distance est d6termin6 en 
955 conf ornament & la relation E87, et 1'on peut poursuivre 

35 en incr<§mentant k (959) jusqu'a la valeur maximum kx (957), 
en fait jusqu'& l'epuisement de tous les symboles bruit6s 
a(k) obtenus en 800. 

0/ 
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Les traitements d' estimation des parametres inconnus du 
modele El du signal, dans l'hypothese "modulation numerique 
lineaire", sont termines. II reste maintenant a verifier que 
l'hypothese de base est correcte, autrement dit a valider - 
ou non - l'hypothese que l'on est en presence d'une modula- 
tion numerique lineaire (stade 990). 

Les traitements precedents ont permis d'estimer les elements 
rappeles en E91. On calcule maintenant la quantite defxnie 
par 1' expression E92, et l'on compare cette quantite a un 
seuil correspondant a la puissance admissible du bruit. Cette 
puissance admissible correspond a la puissance du bruit o* 
calculee dans la partie "calcul de h pm (520), multiple par 
le rapport entre la frequence d'echantillonnage du signal de 
depart et le rythme de la projection (q/T), et ponderee par 
un coefficient /3 (0 < < D pour tenir compte des imperfec- 
tions sur les estimations precedentes des parametres. La 
valeur de /S est fixee par simulation. 

En bref , quand le signal recu correspond ef f ectivement a une 
modulation numerique lineaire, la difference entre le signal 
regu et la modulation numerique lineaire estimee doit etre un 
signal correspondant a du bruit, d'une puissance faible ,• par 
centre, quand le signal recu n'est pas une modulation 
numerique lineaire, ladite difference n'est alors plus un 
signal de puissance faible pouvant etre vu comme du bruit. 

Les traitements decrits peuvent etre mis en oeuvre sur un 
poste de travail du type PC Pentium, avec des outils logx- 
ciels propre a des traitements de calcul, comme le logical 
MATLAB. Les programmes prepares sous MATLAB peuvent Stre 
utilises tels quels, ou mieux transformes en un langage de 
programmation plus rapide a 1' execution, comme le langage C. 
Des utilitaires de conversion de MATLAB en langage C sont 
disponibles par exemple chez les distributeurs du logiciel 
MATLAB. 

Dans certains cas ou il existe plusieurs signaux dans la 
bande de frequence analysee, determinates par analyse 



spectrale a priori, il peut etre necessaire de traiter le 
signal regu par parties, h chacune desquelles on affectera 
une estimee propre de la frequence porteuse. 

La base de fonctions "rectangulaires" decrite plus haut n'est 
pas limitative. On peut recourir h d'autres jeux de fonctions 
qui satisfassent la condition de definir une base, non 
n«§cessairement orthonormee, h condition de tenir compte de la 
correlation induite sur les coefficients de projection. 

II demeure particuliSrement interessant d'utiliser une base 
dont les fonctions se deduisent les unes des autres par 
decalage temporel. Cela autorise en effet une realisation 
avantageuse de la projection sous forme de filtrage, suivi 
d'un echantillonnage. Cette possibility vaut pour d'autres 
types de base de fonctions: 

- si 1'on prend par exemple pour <p x une fonction en cosinus 
sur<§lev<§ de largeur spectrale 2/T, (les autres fonctions <p L 
en etant deduites comme pr6c6demment ) , les spectres sont 
alors ceux de la figure 10B. L'homme du metier comprend que, 
du fait de la forme rectangulaire du spectre de <p lr on perd 
encore moins de signal par la projection ; 

- on peut aussi prendre la base de fonctions des relations 
E34 et E35. Ceci revient (Figure IOC) & filtrer le signal par 
la nouvelle fonction <p 1 (t) et a echantillonner tous les T/q 
(echantillonnage & q points par symboles, avec q > 2). Comme 
le montre la figure 11, comparativement a la figure 8C, 1' ap- 
proximation de h m (t) par sa projection h pm (t) dans cette 
« base » est meilleure. 

II a ete propose plus haut un traitement par matrices de 
covariance entrelacees, applique & la base de fonctions des 
relations E25 et E26. Ce traitement demeure applicable avec 
toute base de fonctions pour laquelle la projection peut etre 
interpretee comme un filtrage suivi d'un echantillonnage tous 
les T/2 (c'est-a-dire a 2 points par symboles). 

Par contre, lorsque la projection correspond & un filtrage 
suivi d'un echantillonnage a plus de deux points par symbo- 

4 



les, la Demanderesse prevoit une autre version dudit traxte- 
ment, que l'on decrira maintenant, en notant q le nombre de 
points par symboles, avec q > 2 (ou q = 2). 

Apres projection, les nouvelles mesures z(k) (coefficients de 
projection du signal) se modelisent sous la forme des 
relations E50. On definit les vecteurs Z L1 (k) selon la 
relation E51 (comme E41, mais avec q a la place de 2 ) . 

La matrice M Ll peut etre formee et traitee comme precedemment. 
Elle est donnee par 1' expression E52. Cette matrxce M L1 
possede pour les memes raisons que precedemment p.L x valeurs 
propres dont L + L, - 1 sont non nulles. La repartition des 
valeurs propres est celle de la figure 13. La cassure du 
diagramme permet ici aussi d'estimer la longueur L (en nombre 
de T) de la fonctioii h pm (t). 

C-est ensuite qu' intervient la difference : si, comme 
precedemment, connaissant L, on calculait la matrice M L , xl 
severer ait que celle-ci, au lieu d'une seule valeur propre 
nulle, aurait q.L - (2.L - 1) valeurs propres nulles (p.L va- 
leurs propres en tout, dont 2L - 1 non nulles, suivant le 
raisonnement donne plus haut). Quand q = 2 , on a 2L - ( 2L ) 
= 1 valeur propre nulle. De 1' unique vecteur propre assocxS 
a cette valeur propre, on en deduit h pm (t). Quand q > 2, xl 
existe plusieurs vecteurs propres associes aux valeurs 
propres nulles. II n'est alors plus possible d'en deduxre 
h (t) d'une manifere simple. 

Pour estimer les coefficients de h pm (t), la Demanderesse pre- 
voit alors de calculer et d'utiliser les vecteurs propres V if 
de v T1 + t * V a Ll , associes aux valeurs propres nulles de la 
matrice M L1 . Les q.L, composantes d'un tel vecteur (d'xndice 
i) sont notees de la facon indiquee en E55, avec i e {L, + L, 
i , p. Li}. 

On forme alors les matrices P ± definies en E56. II s'agit de 
matrices rectangulaires de dimension q.L * qd^ + L>-1, 
formees de la meme facon que la matrice H permettant de 

0( 
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generer le filtre inverse a h pm (t). Chaque colonne d'une telle 
matrice possede q.L composantes . La premiere colonne est 
obtenue en prenant les q dernieres composantes du vecteur V ± 
considere, et en compliant avec des z6ros. Pour le reste, 
chaque colonne se d^duit de la prec^dente par un decalage de 
q cases vers le bas, en compietant en haut par les composan- 
tes de V ± dans l'ordre, puis par des z#ros lorsque ces 
composantes sont epuisees. Ainsi, a la seconde colonne, on 
prend q composantes de plus vers le haut, et ainsi de suite 
jusqu'a la colonne remplie completement avec des z6ros et les 
q premieres composantes de V ± . Le nombre de lignes des 
matrices P ± <§tant fix6 a q.L, on obtient LI + L - 1 colonnes. 

On forme ensuite la matrice G d^finie par 1' expression E57. 
Elle peut §tre consid6r6e comme une forme approch^e d'une 
matrice dite "de projection dans l'espace bruit". 

On recherche les valeurs propres de la matrice G, puis le 
vecteur propre associe a la valeur propre la plus faible. Les 
composantes de ce vecteur propre sont, dans l'ordre, les 
composantes recherchees de h pm (t) dans la base de fonctions 
consid^ree . 

En effet, la Demanderesse a montr6 que cette matrice a une 
seule valeur propre nulle (hors bruit), et que le vecteur 
propre associ<§ a cette valeur propre obeit a la relation E58. 

La puissance a 2 du bruit est dans ce cas estimee en prenant 
la moyenne des valeurs propres presque nulles, divisee par 
LI. L' estimation du train de symboles c m (k) a partir de z p (t) 
et h pm (t) se fait, elle, de la meme maniere avec q > 2, 
qu'avec q = 2, en construisant la matrice H conformement a 
1' expression E63 en prenant la valeur q utilised. 

Ce proc<§de est pr6f6rentiel dans le cas ou l'on a plus de 2 
points par symboles (q > 2). 11 peut aussi etre appliqu<§ dans 
le cas ou q = 2. On estime dans ce cas h pm (t) directement a 
partir de M L1 sans passer par M L . 



L ' invention 
d§crits . 



n'est pas limitee aux modes de realisation 



Pour certaines applications, il peut etre avantageux d'af fi- 
ner 1' estimation approchee f a de la frequence porteuse f 0 . 
Differentes techniques sont disponibles a cet effet. On peut 
par exemple estimer f a en recherchant la frequence, qui apres 
demodulation du signal par celle-ci, minimise la bande 
quadratique du signal. Si l'on parvient a rendre negligeable 
1'ecart residual, alors le stade d' estimation de Af 0 et la 
demodulation qui lui fait suite deviennent sans ob 3 et. 

On notera que le signal de depart capte et enregistre n'est 
pas tou jours assimilable au signal ef f ectivement emis, dont 
il peut differer pour diverses raisons, par exemple des 
trajets multiples de propagation. Il peut en rSsulter des 
amenagements de certains stades du traitement. Les tracts 
multiples, que l'on peut voir de facon imagee comme des 
echos, modifient la forme de la fonction g(t) et plus 
particulierement augmentent sa longueur. La modification de 
la forme ne change pas 1'algorithme puisque celui-ci estxme 
cette forme sans a priori. Il faut par centre augmenter lea 
bornes super ieures LI et N des longueurs de h pm et c^. Ces 
grandeurs parametrables dans 1'algorithme sont a fixer selon 
la longueur maximale des <§chos supposes. 

L ' invention a ete definie essentiellement en reference a la 
modulation numerique lineaire. Bile peut s'etendre a d'autres 
types de modulation auxquels lea stades de traxtement 
proposes conviendraient , au moins en partie. C'est notamment 
le cas pour des types de modulation apparentee, comme la 
modulation numerique de frequence a phase continue, entre 
autres codages par saut de frequence (ou FSK, pour "Frequency 
Shift Keying"). L'homme du metier comprendra que, dans ce 
> cas, la mod«51isation de la frequence instantanSe repond a une 
expression proche de la formule E10 en remplacant Af 0 par 
zero. Tous lea parametres dans E10 sont dans ce cas des 
nombres r6els. 
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ANNEXES 



rCtHRelSaWgO-kTy 2 *'' 



= Re[z(t)] (ED 



p(t)-»p(t)e w 

e j2irf,i (E3) 
x(t) = z(t)z*(t-x) (xpetit) (E4) 

f s 4f-Af„f +Af,] 

z(t)=X;a(k)g(t-kT)e i2 " Ar «« (E10) 

k 

z(t)=Y.(«U* i2 *"°") (g(t-kT)e j2 ^' kT) ) (EH) 

k 

gm(t) = g(t)ej 2 ' tAf «' t (E12) 
a m (k) = a(k)ej 2ltM okT ( E13) 

z (t)=E a »( k )g»( t - kT ) ^ 

k 

Mk), gm (t)) <* 15 > 

Xc.tkJh.O-H)^ 11 -^)*-^-") ^ 

k k 

*(t) = Zc.(k)h.(t-lrf) (E20> 

k 

g m (t)=i> i u t - iT ) w 



i=0 



*(0 = £c m (k)h m (t-kT) 

k (£22) (E23) 

0.(^ = 20, a m (k-i) 

1=1 
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q,,(t) =rect(t) 

H 

Vie{l,..b} *,« = ^1(1) =9,(1-0-1)1) 

Z P « =Sc m (k)h pm (t-W) 
k 

c ro (k) = ta i a m (k-i) 

Qi ( t ) = rect(t) qeIN\{0,l} 



f z(2k+l) = C m (k)h 1+ C in (k-l)h3 + ... + C m (k-L + l)h 21 .. l 
Vk \z(2k+2) = c.(k)h a+ C.,(k-l)h 4+ ... + C ii (k-L + l)h aI . 



r z(2k + l) 
z(2k + 2) 



M k ) = 



z(2k + 2L, - 1) 
z(2k + 2L,) 



(2.L.X2J-,) k 



(E25) (E26) 



(E30) 



(E34) (E35) 



(E36) 



(E37) 



(E38) 



(E41) 



(E42) 
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M L =2X(k).Z- L (k) 

(2.L»2.L) k 



z(2k + l) 
z(2k + 2) 



Z L (k) = 



z(2k + 2L-l) 
z(2k + 2L) 



V = 



+ »»2L 
" h 2L-l 



-hj 



(E43) 



(E44) 



(E45) 



Vkz Lik (t)= ^(k-L+O.Z.O) 
+ C m( k - L+2 > Z 2 W 



+ c m (k) . Z L (t) 

+ 



+ c 



m 



( k+L l) 



(E47) 
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z,(t) 
z 2 (0 



h (t-(-L+l)T) . rect(t) 
pmV V ' ' [0,L,T[ 

h offl (t-(-L+2)T) . rect(t) 



z L (0 



rect(t) 
[O.L,T[ 



(E48) 



7 M = h^ft-Lj) . rect(t) 

Z L+Ll -,W p»V / [0,L,T[ 



(qk + D = C m (k)h, + C n (k-l)h q+1 + ....C ra (k-L+l)h ( ,. l)L+1 
:(qk + 2) = C m (k)h 2 + C m (k-l)h q+2 +....C m (k-L+l)h (q _ l)L+2 



z(qk 



+ q) = C m (k)h q + C m (k-l)h q+q +...-C m (k-L + l)h qL 



(E50) 



%u 00= 



z(qk + 1) 



z(qk + qL,)_ 



M Li =Z Z L .( k ) Z ^ k > 
(q.L, xq.L,) k 



(ESI) 



(E52) 



(q.L.xi) 



ie {Li + L , qLi) 



(ESS) 
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(q.Lxq.fL.+L)-!) 



^(qL.-q + l)" 



v i (qL 1 -2q + l)" 



0 



'i(qL,) 



(E56) 



Vi(qL.) 
0 



_Vi(qL, -q) 



'.<q)_ 



g=£p,p; 



(E57) 



h 2 



pm 

(q.Lxl) 



(E58) 
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H = 
(q.Lx(2L-l)) 



h((L-l).q + l) h((L-2).q + l) 



h(L.q) 
0 



h((L-l).q) 



h(L.q) 



(h\h) _, h- 

[(2L-I)xql.] 

E(a(k)a(k + p)') « £ a < k > a < k + P>* = 



1 pour p = 0 
0 ailleurs 



N 

P(Z)= X a i Z ' 
i = 0 

P 2 (Z) = P(Z).p-(Z) = fia i Z i lfta;: 
9(Z) = 'Z fZ[c.(kK(k+r)j)z-' 

r = -N V K J 

p(z) = K,nM i z i ) 



p 2 (z)=K I .K;n(i-p i z i )(i-p;z i ) 

i = l 



v=^sLwi-isSk(k)l 

k L k 



h(i) 



Kq). 



(E63) 



(E64) 



(E70) 



(E71) 



(E72) 



(E73) 



(E75) 



(E76) 



(E80) 
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E=(a(k)a(k + p)') = 
avec e entier pas trop grand 



X quelconquepour -e<p<e (e e IN et petit) 



0 ailleurs 



(E81) 



1-11 -!^ 2 



a est une constante dont on a fix£ la valeur par simulation 



(E82) 



1 ^, 1 1 



2o* 



(E82) 



Jd 2 (x,y)dxdy 



Jd n (x, y)dxdy <-iy( Jd(x,y)dxdy- Jd„ (x,y)dxdy) A 



(ESS) 



(E86) 



e(k) = eT(||a 



M- e 'f) • 

e ; e{e,,e 2t ...| ensemble estim£ en 900 



3 



(E87) 
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z(t)=£e k g(t-kT)e^ 

k 

(e k ) est le train de symboles estimes 



— est le rythme estim6 
T 

h rat la fonctionde support le plus petit, esttmfc 

pm 

(a,) est le residu de filtrage estim6 
f - f + M n estlaporteuse estimee 



72° 



(E91) 



(E92) 
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Revendications 

1. Dispositif d'aide & 1' analyse de signaux contenant une 
modulation numerique par des symboles, ce dispositif compre- 

5 nant : 

- une mSmoire de signal (100), pour stocker un signal 
numerique complexe (z(t)), repr<§sentatif en amplitude et phase 
d'un signal capte, sur une duree choisie, et 

- des moyens de traitement (200-900) f agences pour rechercher 
10 dans le signal complexe des proprietes relatives k sa 

frequence porteuse et & la modulation de celle-ci, en 
particulier en fonction d'un module de modulation choisi, 
caracterise en ce que les moyens de traitement comprennent: 

- des moyens (200) pour determiner une estimee du rythme (1/T) 
15 de la modulation, 

- des moyens de projection (510), agences pour calculer les 
composantes (z p (t)) du signal complexe dans une base de 
fonctions ( 501 ; cp ± ( t ) ) , laquelle est param6tr6e selon ledit 
rythme (1/T) de la modulation, et 

20 - des moyens de calcul (520), operant sur ces composantes, 
afin de determiner au moins une estimee relative h au moins 
une propriete du signal complexe, dans le groupe de proprietes 
comprenant la forme d" impulsion eiementaire (g(t)) du signal 
complexe, la suite des symboles (a(k)) du signal complexe et 

25 la porteuse (f Q ) dudit signal complexe, 

2. Dispositif selon la revendication 1, caracterise en ce 
que ladite base de fonctions (501) est parametree pour 
presenter au moins deux 6chantillons par periode (T) de la 

30 modulation (q^2). 

3. Dispositif selon la revendication 2, caracterise en ce 
que ladite base de fonctions (501) comprend au moins deux 
fonctions (cp x (t) ,cp b (t ) ) se deduisant I'une de 1 ■ autre par une 

35 translation temporelle de periode choisie (T/2;T/q, q^2). 

4. Dispositif selon l'une des revendications 1 h 3, 
caracterise en ce que ladite base de fonctions (501) comprend 
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des fonctions rectangulaires tempore llement adjacentes les 
unes aux autres . 

5. Dispositif selon l'une des revendications 1 h 3, 
5 caracterise en ce que ladite base de fonctions (501) comprend 
des fonctions du type "cosinus sureieve" • 

6o Dispositif selon l'une des revendications 1 St 5, 
caracterise en ce que les moyens de projection (520) compren- 
10 nents 

- des moyens (521) definissant un filtre numerique, ayant une 
reponse impulsionnelle sensiblement egale k l'une desdites 
fonctions, ce filtre numerique recevant le signal complexe, 
et 

15 - des moyens (523) d' echantillonnage numerique r6p6te de la 
sortie de ce filtre k une cadence choisie (2/T;q/T)o 

7. Dispositif selon l'une des revendications 1 h 6, 
caracterise en ce que les moyens de traitement sont agences 

20 pour determiner (300) une estimee approchee f a de la porteuse 
f Q dudit signal complexe, ainsi que pour demoduler (400) ce 
signal complexe par cette estimee f a , et en ce que les moyens 
de projection (510) sont agences pour op6rer sur le signal 
complexe ( z(t) ) aprfes sa demodulation par cette estimee 

25 approchee, tandis que ladite base de fonctions (501) est de 
basse frequence, comme le spectre du signal complexe apr£s 
demodulation „ 

8. Dispositif selon l'une des revendications 1 & 7, caracte- 
30 rise en ce que les moyens de calcul (520) comprennent des 

moyens de calcul matriciel sur lesdites composantes. 

9. Dispositif selon l'une des revendications 1 & 8, caracte- 
rise en ce que les moyens de calcul (520) sont agences pour 

35 calculer une estimee de la forme d • impulsion eiementaire 
moduiee (9 m ( t )) sous la forme d'une fonction de support 
minimal. 
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10. Dispositif selon la revendication 9, caracterise en ce 
que les moyens de calcul (520) d'une estimee de la forme 
d' impulsion eiementaire moduiee sont agenc6s pour rechercher 
celle-ci comme une fonction de support mimimal, comprenant une 
forme d' impulsion de support minimal ( h pm (' t ))f et un train de 
symboles (c m (k)) associe h ladite forme d' impulsion de support 
minimal . 

11. Dispositif selon la revendication 10 , caracterise en ce 
que les moyens de calcul (520) sont agences pour determiner 
une representation d'un sous-espace de dimension minimale de 
l'espace de fonctions, sous-espace qui contient ledit signal 
complexe, et pour rechercher dans ce sous-espace une direction 
orthogonale Si chacune des tranches du signal complexe, les 
composantes du vecteur propre (V) associe & cette direction 
etant representatives de ladite forme d' impulsion de support 
minimal . 

12. Dispositif selon l'une des revendications 10 et 11, 
caracterise en ce que les moyens de calcul (520) sont agences 
pour determiner en outre le train de symboles (c m (k)) associe 
Si ladite forme d' impulsion de support minimal. 

13. Dispositif selon l'une des revendications 11 et 12 f 
caracterise en ce que les moyens de calcul comprennent : 

- des moyens (610) pour engendrer un filtrage inverse de 
ladite forme d' impulsion de support minimal, puis pour 
appliquer ce filtrage inverse aux composantes du signal 
complexe dans ladite base de fonctions, et 

- des moyens (650) de resolution polyn6miale sur le r6sultat 
de ce filtrage inverse, avec des moyens (659) pour seiection- 
ner un jeu de solutions (^(k)) repondant & des contraintes 
choisies. 

14. Dispositif selon la revendication 13, caracterise en ce 
que les contraintes choisies (659) sont la decorreiation des 
symboles et la minimisation de la variance du module des 
symboles . 



15. Dispositif selon l'une des revendications 10 a 14, 
caracteris6 en ce que les moyens de calcul sont agences (700) 
pour determiner une correction (Af Q ) de 1 • estimation initiale 
(f ) de la porteuse, par recherche d'une frequence de demodu- 
lation du r^sidu qui conduise a une density de probability a 
entropie minimale . 

16. Dispositif selon la revendication 15, caracteris6 en ce 
que les moyens de calcul sont agences (800) pour produire une 
representation du signal complexe demodulS par 1' estimation 
de porteuse corrigee, avec la forme d' impulsion de support 
mimimal . 

17. Dispositif selon la revendication 16, caracterise en ce 
qu'il comprend des moyens (900,950) de calcul d'une estimee 
de 1- ensemble des <§tats possibles (e ± ) des symboles dans le 
signal complexe. 

18. Dispositif selon la revendication 17, caracterisS en ce 
que les moyens (900) de calcul d'une estimee de 1' ensemble des 
Stats possibles des symboles dans le signal complexe operent 
par recherche de maxima locaux de la density de probability 
des symboles. 

19. Dispositif selon la revendication 18, caracteris6 en ce 
qu'il comprend des moyens (950) de calcul d'une estimee du 
train de symboles (a(k)) ef f ectivement present dans le signal 
complexe, par recherche, pour chaque symbole pris individuel- 
lement, de l'6tat possible le plus proche. 

20. Dispositif selon la revendication 19, caracteris<§ en ce 
qu'il comprend des moyens (990) pour reconstruire localement 
un signal ayant la porteuse, la modulation, et les symboles 
estimes, et pour comparer ce signal local au signal initial, 
tel que memorise . 

21. Dispositif selon l'une des revendications prec^dentes, 
caracterise en ce que la^tease^de fonctions (501) est associee 
a 1' ensemble des modul^tijzfn); njumeriques lineaires 
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